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r^ Abstract. In terms of the fully relativistic screened Korringa-Kohn-Rostoker method 

^ we investigate the effect of stacking faults on the magnetic properties of hexagonal 

I close-packed cobalt. In particular, we consider the formation energy and the effect 

V^ on the magnetocrystalline anisotropy energy (MAE) of four different stacking faults 

C^ in hep cobalt - an intrinsic growth fault, an intrinsic deformation fault, an extrinsic 

• fault and a twin-like fault. We find that the intrinsic growth fault has the lowest 

Cd formation energy, in good agreement with previous first-principles calculations. With 

Q the exception of the intrinsic deformation fault which has a positive impact on the 

MAE, we find that the presence of a stacking fault generally reduces the MAE of bulk 

r-f Co. Finally, we consider a pair of intrinsic growth faults and find that their effect on 

O the MAE is not additive, but synergic. 
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1. Introduction 

Within the magnetic recording industry, cobalt alloys such as CoPt and CoPd are of 
great interest due to their large magnetocrystalline anisotropy energies (MAE) [1]. For 
the purpose of magnetic recording, a large MAE of the recording medium is crucial in 
order to maintain stability of the written information as larger areal information stor- 
age densities require smaller grain sizes [2]. In close-packed metals and alloys, stacking 
faults are known to form relatively easily [3] . This is one of the contributing factors to 
the relatively large ductility and malleability that are observed in many such materi- 
als [5]. For a magnetic recording medium, the presence of stacking faults is generally 
considered detrimental, as disturbances in the microstructure will generally worsen the 
signal-to- noise ratio of the medium [3] . Stacking faults may also break the local lattice 
symmetry and, therefore, drastically influence the MAE. 

Experimentally, the effects of stacking faults are generally measured in terms of the 
stacking fault density, which can be determined from X-ray diffraction spectra (see 
e.g. [H SI E])- There are a large number of experimental studies into stacking fault 
formation energies and the effect of the stacking fault density on the magnetocrystalline 
anisotropy for various magnetic recording alloys [H |6]. However, in experiment, the 
real effect of a stacking fault might be obscured by other phenomena, such as migration 
of impurities along the stacking fault, synergies of closely spaced stacking faults, etc. 
Consequently, a number of theoretical methods have been developed for determining 
the properties and effects of stacking faults, see e.g. [7]. In particular, there is a large 
number of first-principles studies of the formation energies of given types of stacking 
faults in metals [Sj [8]. It has been suggested that stacking fault formation energies 
determined from first-principles may be more accurate than experimental measurements 
[3] as theoretical calculations separate the formation energy from any other correlated 
effects on the total energy. The effect on the MAE of a particular stacking fault is, 
however, less commonly explored. In this work, we aim to determine from first principles 
the effect on the MAE of four different types of stacking faults in hep cobalt. 

2. The stacking faults 

Hexagonal planes can be packed either in an ...ABAB... sequence, yielding a hexagonal 
close-packed lattice structure, or in an ...ABCABC... sequence, yielding a face-centred 
cubic lattice structure |9]. In the hexagonal close-packed lattice structure, the stacking 
direction corresponds to the (0001) axis of the lattice, whereas for the face-centred cubic 
lattice structure, the stacking direction is parallel to the (111) axis of the lattice. In a 
hep lattice, a stacking fault is defined as an interruption in the ...ABAB... stacking of the 
hexagonal planes. While there are of course any number of conceivable stacking faults, 
their varying degrees of formation energies and formation mechanisms mean they have 
different probabilities of occurrence |H]. In line with previous work [31 HH], we consider 
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the following four different stacking faults, denoted in standard notation as Ii, I2, E and 
T2 [ni[12]: 

• Ii (intrinsic): •■•BABABCBCB--- 

• I2 (intrinsic): ■••ABAB|CACA--- 

• E (extrinsic): ■••ABABCABAB--- 

• T2 (twin-like): ■■■ABABCBABA--- 

Here the bold face letters or vertical line denote the plane of reflection symmetry of the 
stacking fault. In an intrinsie stacking fault (Ii and I2), the stacking fault is simply a 
shift of one lattice parameter and the stacking on either side is correct all the way up 
to the very fault [9]. The stacking fault Ii is a growth fault while the stacking fault I2 
is a deformation fault [3]. In the extrinsie stacking fault (E), a plane has been inserted 
so that it is incorrectly stacked with respect to the planes on either side of it O US] • 
In a twinASke fault (T2), the stacking sequence is reflected in the fault layer [3]. In the 
following, we refer to the centre of reflection symmetry as the zeroth layer. The two 
layers adjacent to the centre of reflection symmetry are then indexed ±1, and so on. 
Note that in the case of a stacking fault of type I2, the plane of reflection symmetry lies 
in between two atomic layers. Therefore, in the following, for type I2 the atomic layers 
will be labelled by ±|, ±|, . . ., rather than by 0, ±1, ±2, ... as for the other types of 
stacking faults. 

3. Computational details 

For our theoretical study we employed the fully relativistic Screened Korringa-Kohn- 
Rostoker (SKKR) method, in which the Kohn-Sham scheme is performed in terms of 
the Green's function of the system (rather than the wavef unctions) and the treatment 
of extended layered systems is particularly efficient [HI [151 [E] • We used the local spin 
density approximation (LSDA) of density functional theory (DFT) as parametrised by 
Vosko and co-workers [I7j. The effective potentials and fields were treated within the 
atomic sphere approximation (ASA) and an angular momentum cut-off of d-max = 2 was 
used. The magnetocrystalline anisotropy energy (MAE) is calculated using the mag- 
netic force theorem [18] , within which the MAE is defined as the difference in the band 
energy of the system when magnetised along the easy axis (0001) and perpendicular to 
the easy axis. Alternatively, the uniaxial MAE can be calculated from the derivative 
of the band energy, for more details see Ref . [IH] . Only when calculating the MAE, we 
used an angular momentum cut-off of Imax = 3. 

The LSDA+ASA fails in describing the orbital moment and the MAE of Co accurately. 
Similar to our previous work [19] we, therefore, employed the orbital polarisation (OP) 
correction [2011211122], as implemented within the KKR method by Ebert and Battocletti 
[23] . Note that the OP correction was applied only for the £ = 2 orbitals. Excluding 
the OP correction we obtained an easy-plane magnetisation and a MAE of 6.7 /zeV per 
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cobalt atom, while including the OP correction we obtained an easy axis perpendicular 
to the hexagonal cobalt planes and a MAE of 84.4 /^eV per cobalt atom. This is in good 
agreement with the experimental value of 65.5 /xeV [21] and with the experimental easy 
axis being parallel to the (0001) direction. Our result also compares well with that of 
Trygg et al. [2S], who calculated K = 110 /xeV for hep cobalt using a full-potential 
LMTO method including the OP correction. 

In this study we consider an infinite cobalt system, consisting of two semi-infinite bulk 
cobalt systems and an internal region (region J). Region I contains the stacking fault 
and is positioned in between the semi-infinite regions. The combined system is periodic 
and infinite in the plane normal to the (0001) direction. Due to the long-ranged nature 



of the stacking fault effects on the MAE (see section 4.3), the region / in this study 
had to be kept at a size of around 80 atomic layers. More specifically, for stacking 
faults Ii and I2, systems of 80 atomic layers were required, while for stacking faults 
E and T2, 74 atomic layers were required. In order to keep the calculations tractable 
we limited the self-consistent calculations only for a number of atomic layers near the 
stacking fault, and then appended the bulk potentials for the atomic layers further away 
from the stacking fault. We found that it was sufficient to treat only the 20 centremost 
layers self-consistently. In line with previous first-principles studies of stacking faults in 
close-packed metals, we ignored any atomic and volume relaxations (see e.g. [8]). The 
effects of such relaxations are normally negligible because atoms in the faulted part of 
the system tend to retain their close-packed coordination numbers despite the presence 
of the fault [3l |13l [26l ETJ EH |29j. Throughout, therefore, we have used the experimental 
lattice parameter for cobalt, a = 2.507 A. 

4. Results 

4-1- Stacking Fault Formation Energies 

Before exploring how the stacking faults infiuence the MAE of bulk Co, we would like to 
gain an idea of their formation energy. Within the SKKR-ASA scheme, the LSDA total 
energy can be cast into contributions related to individual atomic cells, Ei, comprising 
the kinetic energy, the intracell Hartree energy and the exchange-correlation energy, 
and into the two-cell Madelung (or intercell Hartree) energy, Euad [S]- For a simple 
bulk metal, like hep Co, ii^Mad is, in practice, negligible, while in the presence of stack- 
ing faults it gives a non-negligible contribution due to charge redistributions. However, 
from our self-consistent calculations we found that -^Mad is in the order of 0.1 — 0.2 meV 
per stacking fault. Since the typical formation energy of stacking faults are by about 
two orders higher in magnitude than this value, in the following we consider only the 
layer-resolved (cell-like) contributions to the the total energy. In order to check these 
contributions for artefacts of the appending of the bulk potential, we compare Ei for 
i = —10 (layer with effective potential from a self-consistent stacking fault calculation) 
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to Ei for i = 10 (layer with appended bulk potential), since due to the mirror symme- 
try, these two contributions should be identical. Reassuringly enough, they agreed to 
within a relative error of 10~^, which is well within intrinsic and numerical error of our 
computational method. 

The layer-resolved contributions to the total energy across the systems containing the 
stacking faults Ii, I2, E and T2 is shown in Fig. [1} Herein we observe the expected 
mirror symmetry and that the layer-resolved total energy contributions approach the 
bulk total energy, Eco = —37839.459 eV, towards the edges of each system. From this 
figure it is obvious that the deviation of E^ from Eco is significant up to about 15 layers 
away from the centre of stacking fault. 
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Figure 1. The layer-resolved contributions to the total energy in four hep cobalt 
systems, each exhibiting one of the four different types of stacking fault. The label 
refers to the plane of mirror symmetry. Solid lines serve as a guide for the eyes. 
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The stacking fault formation energy is defined as the difference in the total energy caused 
by the presence of the stacking fault. In order to ensure we include enough atomic layers 
on either side of the stacking fault, we consider the cumulative sums: 

N 

AE(i,,E,T,)(iV) = Y,E^-{2N + l)Eco , (1) 



-N 



and 



N-h 



AEi,{N)= J2 E,-2NEco. 

-N+}: 



(2) 



The formation energy of a given stacking fault X = Ii, I2, E, T2, E\J^, is then defined 



as 



EfoL = lim (AExiN)) . 
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Figure 2. The cell-like part of the formation energy, AEx{N), see Eqs. (IT]) and (l2|, 
of stacking faults Ii, I2, E and T2 in hep cobalt, displayed as a function of the number 
of layers, N, considered in the system on either side of the stacking fault. Solid lines 
serve as guides for the eyes. 



The calculated values of AEx{N) are shown in Fig. |2j Quite obviously, for all types of 
stacking faults, nearly 30 atomic layers (i.e., 15 layers on either side of the stacking fault) 
are required in order to obtain well-converged stacking fault formation energies. The 
fact that the effect of the stacking fault is relatively long-ranged could have significant 
impact on nano-sized systems as the formation energy and, consequently, the likelihood 
of occurrence of a stacking fault could be different depending on its location in relation 
to, e.g., other imperfections as well as surfaces or interfaces in the sample. We obtain 
the following formation energies, with a possible error of ~ 0.1 — 0.2 meV due to the 
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Madelung energy not being included: 

E\^o}m ~ 16 meV ^ 40 mJ ■ m"^ 

E^^orm ~ 48 meV ^ 122 mJ ■ m'^ 

^fJrm ~ 62 meV ^ 160 niJ ■ m'^ 

Ef^m ~ 39 meV ^ 100 mJ ■ m'^ . 

As expected, all stacking faults incur a positive change in the total energy. Of the four 
types of stacking faults considered here, the intrinsic stacking fault Ii has the lowest 
formation energy and the stacking fault E exhibits the highest one. While there is no 
available experiment in literature, the overall results agree well with e.g. Ref. [8]: the 
extrinsic stacking fault formation energy for the close-packed fee metals in this study is 
generally significantly larger than that of the intrinsic and twin faults. Moreover, our 
calculated values for the hep Co growth stacking fault Ii and the hep Co extrinsic fault 
E are close to those obtained by Crampin and co-workers for Ni (which is next to Co in 
the periodic table) [8] : 28 mJ ■ m~^ for the intrinsic stacking fault and 180 mJ ■ m~^ for 
the extrinsic fault. 

4-2. Layer- Resolved Contributions to the Magnetocrystalline Anisotropy Energy 

Because it is calculated directly from the band energy, the MAE can naturally be re- 
solved into layer-dependent contributions, Dj, see Ref. [12]. These layer-resolved con- 
tributions are depicted in Fig. [3] for the different types of stacking faults. Note that 
the mirror symmetry is well reproduced in the layer-resolved MAE contributions for 
all stacking faults. Moreover, the MAE approaches the bulk MAE, Kco = 84.4 /leV, 
towards the edges of all four systems. For stacking faults of type Ii, I2 and T2, the MAE 
contributions become negative at the centre of the fault, favoring thus an in-plane easy 
axis in these layers. This could indicate that these types of stacking faults may act as 
pinning sites. For the type E stacking fault, the layer-resolved MAE contributions near 
the centre are also reduced, retaining, however, very small positive values. 

Furthermore, we note that all stacking faults induce long-ranged oscillations in the 
MAE. For layers of about \i\ > 15, the four stacking faults exhibit very similar trends 
in the layer-resolved MAE contributions. In other words, at about 15 layers away from 
the stacking fault, the presence of a stacking fault still influences the MAE, while the 
particular type of the stacking fault is less significant. This will, however, obviously 
depend on the size of the sample. 

4-3. Finite Size Effect on the Magnetocrystalline Anisotropy Energy 

The long-ranged oscillations in the MAE could cause significant finite-size effects in the 
experimental determination of the MAE of nano-sized samples. We therefore consider 
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Figure 3. Calculated layer-resolved contributions to the MAE for stacking faults Ii 
and I2 (upper panel) and E and T2 (lower panel). The horizontal line refers to the 
bulk MAE, 84.4 /i eV/atom. Solid lines serve as guide for the eyes. 



the following cumulative sums, 
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i^(i„E,T,)(iV) = 5^ A - {2N+l)Kco , 
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and 
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Ki,iN)= J2 D,-2NKco 
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where the MAE of the stacking fault systems of finite width is related to the MAE of 
hep Co. 
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Figure 4. Cumulative sums of layer-resolved contributions to the MAE, Kx(N) {X 
Ii, I2, E and T2), see Eqs. Q and ^. Solid lines serve as guide for the eyes. 



Fig. shows Kx (N) for the four different stacking faults as a function of A^. Surprisingly, 
for N > 5 the I2 type stacking fault appears to increase the MAE, i.e., to strengthen the 
easy axis (0001) {positive effect). As seen from Fig. |3| this is due to the positive MAE 
contributions induced by the stacking fault on neighbouring layers. These apparently 
outweigh the strongly negative MAE contributions induced in the centre of stacking 
fault type I2. This is an unexpected result as stacking faults are typically reported to 
lower the MAE (see e.g. jj|). It should be noted, however, that, of the stacking faults 
studied here, type I2 has the next highest formation energy and it is therefore less likely 
to occur in an equilibrated sample. For stacking faults of types Ii, E and T2, the overall 
change in the MAE with respect to hep Co is negative [negative effect). As noted earlier, 
in the vicinity of these stacking faults, the easy dirction is rotated normal to the (0001) 
axis. This is consistent with the reduction in the total MAE observed experimentally 
by Sokalski et al. in|4|. 



It is quite a remarkable feature that, as seen from Fig. |4| the layer-resolved MAE 
contributions do not settle until at about approximately 35 layers on either side of 
the stacking fault. This long-ranged behaviour could give rise to significant finite-size 
effects in nano-sized samples. Moreover, this might have consequences for theoretical 
investigations into the formation and effects of stacking faults on magnetic properties. 
Typically, in Monte Carlo simulations of stacking faults, interactions between stacking 
faults is kept to around three neighbouring planes Pj. In light of our results, this appears 
to be an uncertain assumption. 
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4.4- Magnetoerystalline Anisosotropy of a Composite Stacking Fault 

Experimentally, the presence of stacking faults is normally quantified in terms of the 
stacking fault density, which is partly a measure of how close the stacking faults are 
located. As the simplest assumption, the change in the MAE due to the presence of 
a number stacking faults in a sample is approximated by the sum of the changes in 
the MAE due to each individual stacking fault. If this were the case, the effect of an 
isolated stacking fault on the MAE could quite straightforwardly be transformed into 
the change in MAE as a function of the stacking fault density. However, the long-ranged 
oscillations in the MAE caused by the presence of each stacking fault indicates that the 
situation is far more complex. 

In particular, we considered two stacking faults of type Ii, separated by three atomic 
layers. In other words, the system exhibits the composite stacking fault: 

■■■ABABCBCBABA--- 

Note that by removing one of the two C-B pairs of atomic layers, a twin-like stacking 
fault T2 is obtained. We have chosen three atomic layers between the centres of the 
two stacking faults, since in dynamical models it is often used as the distance beyond 
which the interaction between stacking faults is neglected (see e.g. Ref. |1]). Moreover, 
we deal with a pair of Ii type stacking faults because this type of stacking fault has the 
lowest formation energy and is, therefore, expected to occur more commonly than the 
other three types of stacking faults. 

The difference between the layer-resolved MAE contributions and the MAE of bulk hep 
Co, 

ADf^'^^ = Df^'^^ - Kco , (6) 

is shown in Fig. [5] for the composite stacking fault. As a comparison, we also show the 
average deviations in the layer-resolved MAE contributions from the bulk MAE of two 
independent type Ii stacking faults, 

ABf.«.> = l(B2^ + D<i.)) -/<■,„. (7) 

If ADl ^ and ADJ ^ were equal for each atomic layer i, the change of the MAE 
due to the presence of the composite stacking fault would be exactly twice that of a 
single Ii stacking fault. However, as shown in Fig. p| AD\ ^ and AD\ ^ deviate 
significantly, particularly in the layers \i\ < 2, i.e., in the layers between the two stacking 
faults. Beyond \i\ > 3, the magnitudes of the MAE contributions are similar for ADJ ^ 
and ADl ^ 5 t)ut with a phase shift of approximately one layer. 



Effect of stacking faults on the magnetocrystalline anisotropy of hep Co: a first-prineiples studyll 



> 
<] 



0.10 



0.05 



0.00 



-0.05 



-0.10 



-0.15 



-0.20 



■ ... J 


. \fm\ . . Mi..^. 


\i 

H If A ndl+Il) 




W T'(i,i,) . 


, 





-40 -30 -20 



-10 10 20 
atomic layer, i 



30 



40 



Figure 5. • : Calculated deviations in the layer-resolved MAE contributions, AD\ ^ ^ , 
of the composite stacking fault from the bulk Co MAE, see Eq. (|6|, and -I- : the 
corresponding average deviations, ADl ^ ^ , of two superposed Ii type stacking faults 
centred on atomic layers i ± 2, see Eq. (l7|. Solid lines serve as guide for the eyes. 



Similar to the case of single stacking faults, we calculate the cumulative sum of the 
MAE contributions of the composite stacking fault, 



N 



Ki,i.{N)= Y.Df^'^'^-{2N + l)Kco 



(8) 



-TV 



and plot it in Fig. [6j Apparently, Ki^i^{N) converges to approximately —1.18 meV 
for large A^, which is almost three times the change of the MAE of the single type 
Ii stacking fault (~ —0.40 meV, see Fig. |4]). Also shown in Fig. |6] is the difference 
Ki^i^{N) — 2Ki-^{N), which appears to settle at approximately —0.38 meV. In other 
words, the two stacking faults interact to yield a stronger negative effect on the total 
MAE as compared to two isolated type Ii stacking faults. This appears to be mainly 
due to MAE contributions from the atomic layers located in between the two type Ii 
stacking faults. This could have significant consequences for predicting the resulting 
MAE in dynamical models used to explain experimental data. To draw any definite 
conclusions, a systematic study of the stacking fault types and separations would be 
required. We expect that such a study would be computationally extremely intensive as 
interlayers (or supercells) of up to approximately 160 atomic layers would be required 
in order to reach the limit in which the two stacking faults are far enough apart not to 
interact. 



5. Summary and Conclusions 

Using the fully relativistic screened Korringa-Kohn-Rostoker method, we have studied 
the MAE of bulk hep cobalt in the vicinity of four different types of stacking faults. 
We find that, in accordance with experiment, most stacking faults have a detrimen- 
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-ft'iiii(^), as defined in Eq. [8|. • : Interaction term of the two stacking faults in the 
MAE, Kj^j^ {N) - 2Kj^ (TV), see. Eqs. Q and ^ Solid hnes serve as guide for the eyes. 



tal overall effect on the MAE. The one exception to this overall conclusion is the type 
I2 intrinsic stacking fault, which, however, exhibits a relatively high formation energy 
and which may, consequently, occur relatively infrequently under standard experimen- 
tal conditions. The effect of a stacking fault on the layer-resolved contributions to the 
MAE is long-ranged, in the order of 15 atomic layers on either side of each stacking 
fault. Motivated by this observation, we investigated a particular composite stacking 
fault and concluded that the MAE of the composite stacking fault is not identical to 
the sum of the MAE of the two isolated stacking faults. A further challenging study 
is proposed regarding the dependence of the 'interaction' of two stacking faults on the 
separation between them. 
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